www.gusucode.com > Matlab在化学工程中的应用 > Matlab在化学工程中的应用/实用化工计算机模拟－Matlab在化学工程中的应用/Examples/Chapter 4/MembraneReactor.m

    function MembraneReactor
% 膜催化反应器的模拟计算
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/06/28 $

clear all
clc
FA0 = 10;
FB0 = 0;
FC0 = 0;
Vspan = [0, 500];
F0 = [FA0, FB0, FC0];
[V,F] = ode45(@ModelEqs,Vspan,F0);
plot(V,F(:,1),'k--',V,F(:,2),'r-',V,F(:,3),'b-.')
xlabel('Volume (liter)')
ylabel('Component flow rate (mol/min)')
legend('C_3H_8','H_2','C_3H_6')

% ------------------------------------------------------------------
function dFdV = ModelEqs(V,F)
Fa = F(1);
Fb = F(2);
Fc = F(3);
k = 0.7;
kc = 0.2;
Cto = 0.2;
Ft = Fa + Fb + Fc;
Kc = 0.05;
ra = -k*Cto*((Fa/Ft) - Cto/Kc*(Fb/Ft)*(Fc/Ft));

dFdV = [ra; -ra-kc*Cto*(Fb/Ft); -ra];